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Abstract. A study of the floating of a circular cone shaped buoy in an ideal fluid 
has revealed some new interesting results. Using reduced variables it is shown, that at 
a crossover value (|) of the ratio of the specific masses of the fluid and of the buoy, 
the anharmonicity of the oscillation is the highest and that, unexpectedly, above this 
crossover value the normalized period is constant. 



1. Introduction 

Anharmonic oscillators are useful examples of nonlinear phenomena. Many vibrating 
systems found in the real world are nonlinear whether they be macroscopic mechanical 
oscillators [U 121 El H] or microscopic atomic oscillators^. The pendulum at high 
angles is a classic example of an anharmonic oscillator [6l U\- A particular feature is 
that the restoring force is equivalent to a spring that softens at large amplitudes i.e. 
F = —ki6 + /c2^^.[l] In this work we do a detailed analysis of an oscillator materialized 
by the bobbing cone; we assume that the movement is a vertical translation and that 
the top of the solid is always emerged and the base is always submerged in an ideal fluid 
for which the restoring force is due to Archimedes' principle. The restoring force being 
of the type F = —ay + by"^ — cy^ (a, b, c positive parameters) means that the oscillator 
does not move symmetrically about the origin notwithstanding being periodic. 



2. Restoring force 

Figure ([1]) schematically represents a floating cone of radius R, height h and specific 
mass p, partially immersed in an ideal fluid of specific mass pf {p < pi). The position 
of the solid is chosen to be the coordinate of the point P which coincides with the 
origin of the reference axis when the buoy is at equilibrium. This origin (point O) is 
the intersection of the axis of the cone at equilibrium with the plane of the free liquid. 
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Figure 1. The floating cone in two positions: a) equilibrium; b) at an instant t, with 
coordinate y. The referential axis is vertical downward. 



The resultant of the forces acting on the sohd is F = pVc g — Pi (Vq + ^V)g where 
Vq is the volume of the immersed cone in static equilibrium {P = O i.e. y = 0), 14 is 
the volume of the cone, Al^ is the difference between the immersed volume V and Vq 
and g is the gravity acceleration. Therefore: 



pV^g = piVog 



Ay = K 

where a = pj pi <\. 
Then, 



h) 



3(1- a) V3 (^|)' + 3(l-af/=^ (|) 
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where Fq = pfgY^. 

Introducing the reduced variables y — y/h and F — F/Fq, the reduced restoring 
force is: 

F^-[f-3{l- ay/^f + 3(1 - a)^/'y] = -Ps{y, a) 

where Psiy, a) is a third degree polynomial with coefficients that depend on the 
parameter a. 

The reduced coordinates of the vertex and of the center of the base are, at 
equilibrium, respectively: yy — —^1 — a and Ijq = 1 — ^1 — a . Therefore, due to the 
restrictions of the movement —IJq <y< —yy. However, another condition should be 
imposed on the value of the initial position t/o (or its reduced form ^q) as a consequence 
of the potential energy barrier, i.e. Fp(|/o) < E-p^—yy) and Ep{yo) < Ep{—yc). 



3. Potential energy 



The potential energy associated with the restoring force F, assuming that Fp = at 
the equilibrium position {y — 0) is 

3(1 3(1 (1)1* 



ry 
Jo 
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= PigVc r [y' - 3(1 - ay/^'t + 3(1 - ay/''y] dy. 
Jo 

In reduced form, Ep = Ep / Eq = Piijj, a), with 
Eq = pighVc = FqH and 

P,{y,a) = ly^iV) dy = ^t-{l- o^f't + ^(1 - af^'t- 
The choice for the initial values {y = yo and zero velocity) is subjected to the relations 

-Vc <yo< -Vy ^ -1 + <yo< - a 

PiiVo) = min {P4{-yc),P4{-yv)} 

These equations express the condition that the oscillation has the maximum energy 
compatible with the fact that the cone is neither completely immersed nor completely 
emerged. 

It is easy to verify that 

p^{-yc) = p^i-i + v^r^) = ^ (1 - «)^^^ + « - ^ 



Then 



P^i-yy) = P.i^T^) = ^ (1 - a)'/'. 

P,{y,) = min 1^ (1 - «)^/^ ^ (1 - a)'/' + « - ^} ; 

for a = 3/4 =^ Pi{~yc) = P^i—Vv)- So, if the initial velocity is zero then the initial 
amplitude y^ should be chosen along with 

a<l^yo = -yc = -l + v^l - a and Erotai = | (1 - ^Y^^ + " - f 
a> l^yQ = -yy = ^/l - a and Erotai = f (1 - a)^^^- 

The other limit y^ of the interval of y can be obtained from the roots of the polynomial 
-^4(^1) ~ Erotai = 0. There are four solutions (two real and two complex conjugate); the 
real are the relevant solutions: [Eqs. (E])] and y^ [Eqs. ([3])] 

1/3 



-44 + 54a + 6v^-27 + 81(1-q)2 + 30q 



-44 + 54a + 6v^-27 + 81(l-a)2 + 30a 



a>i^yi =-(^-1)^ 



a. 



4 

For each a value there is a definite maximum of total energy and a corresponding interval 
of amplitude [a,b] = [—yQ,y-i] for a < 3/4 and —yy] for a > 3/4 given by Eqs. ([2]) 
and ([3]). It means that the procedure to initiate the movement should be as follows 
for a < 3/4 the base of the cone should be raised near the free surface of the liquid 
for a > 3/4 the vertex of the cone should be lowered until near complete immersion 
for a = 3/4 the choice of the base or vertex to initiate the movement is irrelevant since 
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the amplitude of the oscillation, that has a maximum (equal to h) corresponds to a 
displacement between the base and the vertex i.e. y G [— ^c? ~^v]- This can be seen in 
Fig.®. 
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Figure 2. The reduced total energy, E = E/Eq, and the reduced interval, y ~ y/h, of 
the oscillations as a function of a. The maximum of the energy E = 3(1/4)^/^ occurs 
at the crossover a — 3/4. 



The restoring force (reduced value) F = F/Fq = —P^{lj) and the potential energy 
E-p = Ep/Eq = P^ijj) can be represented in the interval of oscillation with a as 
parameter [Fig. (jS]) and Fig. (j4])]. 




Figure 3. The reduced restoring force as a function of the reduced displacement, 
F = f{y) , within the proper interval of oscillation taking a as a parameter. 
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Figure 4. The reduced potential energy Ep — f{y), witliin ttie proper interval of 
oscillation taking a as a parameter. 



4. General discussion of the type of oscillation 



We use the Newton equation to find the position of the point P of the cone as a 
function of time. The adoption of reduced values continues to offer simplification 
in the final equation. Let us define the unit of time to = a/2 (h/g). Therefore 
t = t/to =^ dy/dt = (h/to) dy/dt. The unit for f is fo = ^/gh/2 which points to 
17= v/vq = dy/dH. On the other hand d'^y/dt^ 
motion in terms of the reduced variables is 

^ , d-'y tl F 



(h/tl) d'^y/di , i.e. the equation of 



Finally 



h d^y 



d^y 2- 

— 2. J? 

df a 



d^ 



h pioVc 



2 F 
a Fq' 



d^y 2 ^ 
di « 



0. 



The system does not oscillate symmetrically about the origin since the force (polynomial 
P3) shows odd and even powers. For small displacements ^ <C 1, Pzijj) ~ 3(1 — 
the equation is approximately linear and so the movement is quasi harmonic with 
period Tso = 2Ti^Ja/Q . — a). In other physical situations the oscillator exhibits 

anharmonic displacements and the period will depend on the amplitude (or total energy). 
The solutions y = y(t) and v = v(i) were obtained by computational methods and are 
represented in Figs. ([5]), (j6]), ([7]) and ([8]). 



5. The Period 



The conservation of the energy of the oscillator allows the determination of the velocity 
as a function of displacement: dv/di + {2/a)P3 = =^ vdv = —{2/ Q.)P^dy =^ v = 





Figure 6. Displacement of the cone (reduced values) as a function of the reduced time 
for several values of the parameter a> 3/4. 



±(2/v^) ^Etotai - Pa. 

( y^3{l - a)4/3 + 4a - 3 - y4 + 4(1 _ a)i/3y3 _ g(i _ ^^2/3y2 f 
^~ \ ^3(1 - -y^ + 4(1 - a)V3y3 _ g(i _ c,)2/3y2 ^^>3^ 

The phase space representation is shown in Fig. (Q. From the symmetry of these 
curves relatively to the horizontal axis, we can conclude that the interval of time between 
two zeros of the velocity or two extremes of the displacements is equivalent to half of 
the period of the movement. Since di = dy/v, T/2 = j^di,- Then, 
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Figure 7. Reduced value of the velocity as a function of the reduced time for several 
values of the parameter a < 3/4. 




Figure 8. Reduced value of the velocity as a function of the reduced time for several 
values of the parameter a > 3/4. 
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The limits of integration and the integrand function have been determined above. 
The result is only a function of a. It is possible to get a close solution for the period 
in terms of complete elliptic integrals of the first kindjl]. For a < 3/4, the analytical 
solution is far too complicated to be considered and so numerical methods were used. 
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Figure 9. The phase plane representation {v as a function of y) for some values of a. 



However, the solution for a > 3/4 was easily obtained: 
T = 




^ + 4(1 - ay/^t - 6(1 - 



C- 



where K (^2/(2(1 + Vs))) is a complete elliptic integral of the first kind of argument 
^2/(2(1 + VS)). The numerical value of C is: C = 3.059908075. An interesting result 
that should be outlined is the fact that the normalized value of the period for a > 3/4, 
Tnorm = T/Tso (Tso, pcriod of small oscillations) is a constant, i.e. does not depend 
on a: 

T„ = VMlEl = ^c= 1.192900269. 

^/T^ 277 271 

The period of the cone oscillations at the highest energy as a function of a is represented 
in Fig. ( ITOl) together with the normalized values. 

So far, we have represented the oscillation with the highest energy or maximum 
value of the initial amplitude satisfying conditions ([1]). It is interesting to examine 
closely the dependence of the period with the energy of the oscillator from zero to 
the highest value given by Eqs. ([2]). This is shown in Fig. (ITT]) taking, as usual a as a 
parameter. It is interesting to observe that the numerically obtained values {E, T) fit 
100 % a second degree polynomial. 
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Figure 10. The reduced (left axis) and normalized (right axis) period as a function 
of a, calculated for the highest value of its total energy. The derivative dT/da is also 
shown in the inset. 
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Figure 11. The variation of the period with the energy (reduced values) using a as a 
parameter. The energy varies from to its maximum value (Eq. ©j- 



6. Fourier analysis. 

The Fourier analysis is a very useful tool to investigate the harmonic components of 
a periodic function of time that derives from a non linear equation, as it happens in 
our study From Figs. ([5]) and ([6]) one can conclude that the displacement y = y(t) 
is a periodic symmetric even function with a non null mean value. The Fourier 

oo oo 

decomposition gives for the general case: y{t) = ao + ^ ctn cos 

(^f ) + E &nsin (nf ); 

n=l n=l 




in this case all bn = since y{t) is an even function. The coefficients a„ were calculated 
by numerical methods once the fundamental frequency is known. Only the first four 
coefficients were significant. The Fourier components for the particular case a = | 
that corresponds to the maximum of the highest value of energy are represented in 
Fig. (1121) . It also exhibits the highest coefficients which indicate the highest degree of 
anharmonicity. The reduced values of all calculated coefficients (n = 0, ...,3) are also 
shown in Tabledl 

Table 1. Values of the coeficients a„ {n = 0, . . . , 3) of the Fourier analysis of the 
reduced displacement, y = y(t). 
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7. Conclusions 

Most real oscillators contain anharmonic components. The present study of a floating 
cone movement is by all means an interesting case of an anharmonic oscillator. The fact 
that the restoring force is a polynomial of the third degree with non null coefficients 
(except the independent term) imply an asymmetry of this force and the related 
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potential curve. Comparing this to well known based mass-spring systems results 
that the equivalent spring stiffness is no longer constant. It varies as a sum of two 
contributions that have opposite sign from a specific position. In solid state physics 
we encounter these kind of forces such as the cohesive force containing a short range 
repulsion (hard sphere interaction) and a long range attraction. The former varies with 
the displacement much faster than the latter. 
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